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ABSTRACT 

Using small automated telescopes in Arizona and Hawaii, the HATNet project has detected an 
object transiting one member of the double star system ADS 16402. This system is a pair of GO 
main-sequence stars with age about 3 Gyr at a distance of ~139 pc and projected separation of 
^1550 AU. The transit signal has a period of 4.46529 days and depth of 0.015 mag. From follow-up 
photometry and spectroscopy, we find that the object is a "hot Jupiter" planet with mass about 0.53 
Mj and radius ~1.36 i?j traveling in an orbit with semimajor axis 0.055 AU and inclination about 
85?9, thus transiting the star at impact parameter 0.74 of the stellar radius. Based on a data set 
spanning three years, ephemerides for the transit center are: Tc = 2453984.397 -f Ntr * 4.46529. The 
planet, designated HAT-P-lb, appears to be at least as large in radius, and smaller in mean density, 
than any previously-known planet. 

Subject headings: stars: individual: ADS 16402 AB - planetary systems: individual: HAT-P-lb - 
binaries 



1. INTRODUCTION AND SUMMARY 

The 10 currently known transiting extrasolar planets 
occupy a special place among the nearly 200 known exo- 
planets because their photometric transits yield unam- 
biguous information on their masses and radii. Five 
of these transiting planets were found through spectro- 
scopic and photometric follow-up of candidates discov- 
ered by the OGLE project (t hese are OGLE-TR-lO b. 
56b, 111b, 113b, and 132b; lUdalski et all I2002albll(]. 
[^03, 2004, respectively). For references to the follow- 
up confirmation papers for these objects, see the Ex- 
trasolar Planets Encyclopedia^^. The parent stars of 
these planets are all rather faint, making effective fol- 
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lowup somewhat difficult. The other five^ transiting 
exoplanets orbit nearby, bright st ars: HD 20 9458b - 
iCharbonneau et alJ ll200a) : iHenrv ll2000l). TrES- 

1 - lAlonso et al.l ll20ol. HD 1490 2 6b - ISato"eran 
120051) ■ HD 189733b - IBouchv et alJ M^ . XO-lb - 
iMcCuUough et all l)2006|) . These planets are of special 
interest, as accurate parameter determination as well as 
other types of follow-up are facilitated. For this reason, a 
number of wide-fiel d surveys, using smaU telesc opes, are 
underway (see, e.g. ICharbonneau et ani2006bj) . Indeed, 
two of the 5 known transiting planets around nearby 
stars, TrES-1 and XO-lb, were detected by such surveys. 
(The other three were first detected by radial velocity 
measurements) . 

The HATNet project,!''^ initiated in 2003 by G. A. B, 
is also a wide-field survey that aims for the discovery 
of transiting planets around bright stars. It currently 
comprises 6 small wide-field automated telescopes, each 
of which monitors 8° x 8° of sky, typically containing 
5000 stars bright enough to permit detection of planetary 
transits through the typically 1% photometric dips they 
induce on their parent stars. The instruments are spread 
in a two-station, longitude-distributed network, with four 
telescopes at the F. L. Whipple Observatory (FLWO) in 
Arizona, and two telescopes at the Submillimeter Array 
(SMA) at Hawaii. Technical aspects of HATNet will be 
described in a forthcoming paper, but for instrumenta- 
tion, observations, and data flow, sec Bak os et alJ l|200^ 
l2?iol . Here we report on the first detection by HAT- 
Net of a transiting extrasolar planet. This planet, which 
we hereafter label as HAT-P-lb, orbits a lOth magni- 
tude star, and thus is well-suited for follow-up with large 
ground-based and space telescopes. 

The existence of HAT-P-lb was first inferred from a 

" Discovery of TrES-2 lO'Donovan et al.l2003) was announced 
during the refereeing of this paper. 
www.hatnet.hu 
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Fig. 1. The left panel shows a 70' X 70' (300 pixels wide) stamp of a typical image taken by HAT-5. The stamp width is about 
l/7th of the original frame. ADS 16402 is in the middle (the two components are merged into a bright star). Peaks of saturated stars are 
masked, and shown in white. The right panel is a stamp from an image taken with the FLWO 1.2 m telescope. It is 6.'8 X 6.'8 (600 pixels 
wide), and corresponds to the small box on the left panel. ADS 16402B(marked with the cross) is well separated from ADS 16402A. 



0.6% (6 mmag) deep transit-like dip seen in the /-band 
light curve of the stel lar system ADS 1640 2 ( Aitkcn Dou- 
ble Star Catalogue; lAitken &: Doolittl3 ri932) obtained 
by the combination of the HAT-5 instrument at FLWO 
and HAT-8 at the SMA. The members of this binary 
stellar system are a pair of nearly identical GOV stars, 
ADS 16402A (J « 9) and ADS 16402B (/ « 9.6) sep- 
arated by 11'.'2. The system appears as a single, elon- 
gated source in the HAT images because this separa- 
tion is less than the 14"pixel^^ image resolution. Ini- 
tial follow-up spectroscopy at the FLWO 1.5 m telescope 
showed the absence of radial velocity variations of either 
star at a level of 1 kms~^ or larger, demonstrating that 
the source of the photometric dips could not be due to a 
common false positive, namely an M dwarf star transiting 
in front of one of the stars. Subsequent spectroscopy at 
the Subaru tel escope using N2K proj ect observing time 
ijFischer et all 120041: ISato et aljf2005l|) . and later at the 
W. M. Keck telescope revealed sinusoidal velocity varia- 
tions of ADS 16402B with an amplitude of 60m s~^, thus 
suggesting the existence of a 0.5A'fj planet transiting that 
star. After excluding blend scenarios, the combination of 
the photometric and spectroscopic data leads to the con- 
clusion that this object, HAT-P-lb, is indeed a planet, 
and has a radius comparable to that of HD 209458b, and 
a somewhat lower mean density. In the remainder of this 
paper we discuss all of the above points in more detail. 

2. OBSERVATIONS 

2.1. HATNet Discovery Photometry 

ADS 16402 (also BD+37 4734, a = 22''57™46?8, S = 
-H38°40'28", J2000.0) is a visual double star system 
that also bears the name HJ1832 , as it was discovered 
by John Herschel l)Hersche]l 1 1 83 1() . The system lies in 
HATNet survey field G205 (centered at a = 22''55", 
S = +37°30'). 2059 observations were made of this field 
by the HAT-5 telescope and 901 by HAT-8 between 2003 
October 13 and 2004 January 30. A small stamp taken 
from a typical HATNet image is shown on the left panel 



of Fig. n Exposure times were 5 min at a cadence of 
5.5 min. Light curves were derived by aperture pho- 
tometry for the 6400 stars in G205 bright enough to 
yield photometric precision of better than 2% (reach- 
ing 0.3% in some cases). In deriving the light curves, 
we made use of the Trend F iltering Algorithm (TFA; 
iKovacs. Bakos. fc NovesI l2005() to correct for spurious 
trends in the data. This was particularly important for 
the light curve for ADS 16402, for which variable blend- 
ing distorts the combined shape of the two close-lying 
unresolved components; indeed the shallow transit sig- 
nal in the combined light of the two stars is only barely 
visible in the raw light curve, while it is readily seen after 
application of TFA. 

We then searched all the light curves for characteristic 
transit signals, using the Box- fit ting Least Squares (BLS; 
IKovacs. Zucker. fc Mazehl 12002(1 algorithm, which seeks 
for box-shaped dips in the parameter space of frequency, 
transit duration, and phase of ingress. Candidate transit 
signals with the highest detection significance were then 
examined individually, to isolate those with the best com- 
bination of stellar type (preferably main-sequence stars 
of spectral type mid-F or later), and depth, shape, and 
duration of transit. One of these was ADS 16402, for 
which we identified a prominent boxcar-like periodicity 
with period 4.4656 days and transit depth in the com- 
bined light of the two stars of 6mmag. Portions of five 
transits were observed during the observing interval. The 
upper panel of Fig.Elshows the combined data for these 
five transits, with different symbols for observations from 
HAT-5 (Arizona) and HAT-8 (Hawaii). 

2.2. Rejection-mode Spectroscopic Observations 

Initial follow-up observ ations were ma de with the CfA 
Digital Speedometer fDS: lLathamll99 2?) in order to char- 
acterize the two individual stars comprising ADS 16402, 
as well as to search for evidence that the transit signal 
could be due to the transit of a late M dwarf in front 
of either of the stars, as has been found for many oth- 
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erwise plausible planet transit candidates in the HAT- 
Net database. These observations indicated that both 
sources indeed have values of Teff and logg consistent 
with main-sequence late F or early G dwarfs. The DS 
data revealed no sign of radial velocity variation in ei- 
ther of the two components, with upper limit of detection 
about 1.7kms~^. For a period of ~4.5d this limit corre- 
sponds to a secondary mass of 15Mj, thus ruling out the 
possibility that the transit signal is due to a low-mass 
stellar companion transiting one of the two stars. 

2.3. Follow-up Photometry Observations 
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Fig. 2. — The upper panel shows the HAT-5 (open circles) and 
HAT-8 (filled circles) TFA corrected light curve from FLWO, Ari- 
zona and SMA, Hawaii, respectively. Altogether five transit events 
were observed successfully between the two sites. The data are 
phased at period P=4. 46529 days. Zero phase corresponds to the 
time of mid-transit using the ephemeris discussed in the text. The 
second and third panels show the follow-up observations taken with 
KeplerCam on the FLWO 1.2 m telescope in z-band in 2006 July 
23 and 2006 Aug 31, respectively. Unfortunately only the central 
and egress portions of the transit were observable. The bottom 
panel shows the /-band data of a full-transit taken by the Konkoly 
Observatory 0.6 m Schmidt telescope at 2006 Sep 5. Superimposed 
on the lower three panels are our best fits (§1^1 to the transit light 



After identification of the ADS 16402 system as a can- 
didate for harboring a transiting extrasolar planet, we 
used the TopHAT telescope to search for evidence of a 
transit. TopHAT is a 0.26 m aperture automated tele- 
scope at FLWO, part of th e HATNet, designe d for follow- 
up of transit candidates l|Bakos et al.ll200'^ . With its 
pixel scale of 2.2"pixel~"'^ it is able to resolve the two stars 
in the system. However, data were taken during the com- 
missioning phase of TopHAT with strongly variable fo- 
cus, and furthermore, conditions were non-photometric. 
Nevertheless, we identified in the data a transit ingress 



event associated only with the star ADS 16402B. 

Recent much higher quality observations using Kepler- 
Cam on the FLWO 1.2 m telescope have confirmed this, 
revealing the egress from the transit on 2006 July 23, 
and later, on 2006 Aug 31. A small stamp from an 
FLWO 1.2 m image that clearly shows the two compo- 
nents of ADS 16402AB is exhibited in the right panel 
of Fig. n The FLWO 1.2 m data were reduced using 
our own aperture photometry software (FiHAx/f iphot, 
Pal et al., in preparation). The photometry on individ- 
ual frames was then transformed to a selected reference 
frame using some one hundred stars as comparison in the 
23' X 23' field provided by KeplerCam (we refer to this 
step as smooth magnitude transformation). The result- 
ing light curve for ADS 16402B already showed a smooth 
transit curve with very high S/N, but low amplitude sys- 
tematic variations were present. As a final correction, we 
exploited the fact that ADS 16402B has an ideal com- 
parison star, namely ADS 16402A, with almost identical 
brightness, spectral type, and spatial position. The data 
shown in the second and third panels of Fig. [21 thus cor- 
respond to the differential magnitude of ADS 16402B to 
ADS 16402A after the smooth magnitude transformation 
procedure (shifted back to an arbitrary z « 9.53 out-of- 
transit level). 

Unfortunately both FLWO 1.2 m transit events were 
partial, thus we seeked for further opportunities to ob- 
serve a full transit. Such an event was observed re- 
cently in /-band, on 2006 Sep 5, using the 60/90/180cm 
Schmidt telescope at the Piszkesteto mountain station 
of Konkoly Observatory. Data are shown on the bottom 
panel of Fig. [2 and were reduced in a similar manner to 
the FLWO 1.2 m data. 

The HATNet, TopHAT, FLWO 1.2 m and Konkoly 
Schmidt data together span 1058 days, or 237 periods. 
Combining all the transit timings in these data we find 
improved ephemerides: Tc{HJD) — Eo + Nn- ■ P, where 
Eo = 2453984.397 ± 0.009, and P = 4.46529 ± 0.00009 
(days) , and Ntr refers to the number of transits since the 
recent Konkoly Schmidt observation that is a full transit 
observation, hence it has the best timing measurement. 

2.4. High-precision Spectroscopic Observations 

In parallel to obtaining photometry follow-up data 
with the FLWO 1.2 m telescope, we obtained precise 
Dopplcr measurements of both ADS 16402A and B with 
the High Dispersi on Spectrograp h (HDS) on the 8.2 m 
Subaru telescope ijNoguchil 120021) to confirm the plane- 
tary nature of the HATNet photometric signal. 

The HDS spectral format spans 3500 - 6IOOA over a 
mosaic of two CCDs, and the 0'.'8 slit yielded a resolu- 
tion (A/ A A) of 55000. With only four observations we 
determined that the brighter component, ADS 16402A, 
had low rms velocity scatter while component B exhib- 
ited velocity variations that were well represented by a 
sine curve with an amplitude of about 60m s""'^, a pe- 
riod equal to the photometric period, and a time of con- 
junction consistent with that predicted from the tran- 
sit ephemeris. To improve phase coverage, nine obser- 
vations of ADS 16402B were subsequently obtained at 
W. M. Keck telescope ove r eight nights in 2006 July using 
HIRES l|Vogt et alJll994j) . The spectrometer sht at Keck 
is 0'.'86 yielding a similar spectral resolution of about 
55000 with a spectral format from - 3200 - 88OOA. The 
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TABLE 1 

Summary of stellar parameters for 
ADS 16402AB. 



Parameter ADS 16402A ADS 16402B 

From SME analysis (§ 

reff(K) 6047 ± 56 5975 ±45 

logg 4.13 ±0.10 4.45 ±0.06 

i;sini(kms-l) 7.1 ±0.3 2.2 ±0.2 

[Fe/H](dcx) +0.12 ±0.05 +0.13 ±0.02 



Simultaneous evolutionary track fittingi 13.21 


Mass (Mq) 


1.16 ± 0.11 


1.12 ±0.09 


Radius (-Rq) 


1 no_i_0.14 


1 -itriO.lO 

Q-OZ 

o-i8±8:lI 


logL*{L0) 


0.26 ± 0.15 


Mj (mag) 


3.4 ±0.3 


3.7 ± 0.3 


Age (Gyr) 


3.6 




Z 


0.025 





Other parameters 



I (mag) 9.035 ± 0.05 9.563 ± 0.05 
A/ (mag) 0.53 ± 0.03 

Mean Dist. (pc) 139±?i 
Proj. Sep.(AU) 1550±lfg 



Subaru and Keck data were used to characterize more 
fully the stellar properties of the system (§0), as well as 
to obtain a definitive radial velocity orbit (§EIl. 

3. STELLAR PROPERTIES OF THE ADS 16402 SYSTEM 

3.1. SME analysis of the individual stars 

We have carried out spectral synthesis modeling of the 
iodine-free Keck template spectra fo r both ADS 16402A 
and B using the SME software ijValenti &: Piskunovl 
I1996(1 . plus the wavelength ranges and atomic line data 
described bylValcnti & Fischer (2005). Results for Toff, 
logg, usini and [Fe/H] for both stars are shown in the 
first four lines of Table Q Given the derived values of Teff 
the stars appear to fit the MK classification of GO, rather 
than the F8 classification reported in the literature. 

3.2. Simultaneous evolutionary track fitting 

The two stars are very likely to comprise a physical 
pair, for the followin g reasons: First, they have com- 
mon proper motion (" Halbwachsl Il98(j) . Second, their 
relative apparent magnitudes compared to the abso- 
lute magnitude expected for main-sequence stars with 
Tcff and log g given in Table ^ indicate that they are 
at a common distance. This also implies that their 
true tangential velocities are very similar, and also that 
their real separation is small (projected separation be- 
ing ^^1550 AU). Third, based on the six DS veloci- 
ties for each star acquired over a time interval of 234 
days, they have common radial velocity within the er- 
rors: Kad.A = -3.43 ± 0.32 (rms) ± 0.14(sys)kms-i, 
V,^d,B = -2.94 ± 0.56 (rms) ±0.25(sys)kms-i. Finally, 
they have identical metallicity within uncertainties. The 
physical companionship implies that the two stars are 
coeval, and provides the opportunity to perform simul- 
taneous evolutionary track fitting. 

Furthermore, because of their common distance, the 
difference in their apparent magnitudes yields an accu- 
rate measurement of the ratio of their luminosities. The 
/ magnitudes of ADS 16402A and ADS 16402B respec- 
tively are 9.03 and 9.56. These values were determined by 
cross-correlating the / magnitudes for the Landolt (.1992i) 
standard stars with the J, H and K magnitudes of the 



2MASS Point Source Catalogue ijSkrutskie et al.ll2006(l . 
and by deriving a weighted linear regression between J, 
H, K and /. The linear relation based on ^ 400 stars is 
very well defined, and has an rms of 0.06mag. The / mag- 
nitude difference. A/ — 0.53 ± 0.03, is more accurately 
determined than the individual magnitudes, because the 
stars are nearly identical in spectral type and brightness 
so that systematic errors are significantly canceled. 
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Fig. 3. — Stellar evolution tracks (solid curved lines) for 
ADS 16402A and ADS 16402B corresponding to the best mutual 
fit (filled circles with error bars). The two stars are physically 
bound, are assumed to be coeval, and to have the same abun- 
dance, as determined by the SME analysis. Initial parameters from 
SME are shown with open circles with their respective (overlap- 
ping) error-boxes. The two tracks stop at 3.6 Gyr, so that the line 
joining their endpoints defines the local isochrone. The straight 
dashed lines plot the ZAMS for metallicity Z=0.020 and Z=0.025. 
For clarity the ZAMS for Z=0.030 is not shown in the figure; it 
is displaced by the same amount above the ZAMS for Z=0.025 as 
the latter is above the ZAMS for Z=0.020. Dotted lines illustrate 
lines of equal radius. The position of the Sun and its track to its 
current age are also shown. 

By exploiting the physical association of the two stars, 
their common abundances, age, and their known lumi- 
nosity ratio, we carried out a simultaneous evolutionary 
track fitting. This in turn leads to a much better deter- 
mination of the parameters of ADS 16402B than would 
be the case if it were a single star. 

Initial parameters (Fig.|21 open circles) and their error 
boxes (solid line rectangles in Fig. O were based on the 
atmospheric parameters obtained with the SME analysis 
of the Subaru spectra. While the errors in Toff are taken 
directly from Tabled the absolute luminosities rely on 
the estimated gravity and iteration with the mass de- 
termination, and this procedure is imprecise. However, 
the ratio of the luminosities of the two stars is well con- 
strained, and is used in initializing the isochrone fitting 
to both stars. 

The value of [Fe/H] found from spectral synthesis 
corresponds to the range of metal abundance by mass 
Z=0. 020-0. 030. We have calculated steUar evolution 
tracks of ADS 16402A and ADS 16402B using Z in this 
range (Fig.|31). The straight Hue that would connect the 
centers of the two error boxes is steeper than the theo- 
retical isochrones (the straight line that would connect 
the end tips of the two evolutionary tracks) for any Z in 
the range Z=0. 020-0. 030. All possible solutions are some 
sort of compromise within the observational error boxes, 
imposed by the requirement for equal age, and with the 
obvious additional constraint that each track must begin 
on the ZAMS appropriate to the chosen value of Z. For 
different Z's, different mass ranges sweep the two error 
boxes - smaller ones for Z=0.020 and progressively larger 
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TABLE 2 
Radial Velocities for ADS 16402B. 



HJD - 245o0U0 


RV 


Uncert. 


Observatory 


(days) 


/ 1 \ 

(m s '-) 


(ms '-) 






-2.89 


5.00 


O U iJdi LI 


899.12678 


52.53 


4.77 


Subaru 


900.12008 


-34.51 


4.65 


Subaru 


901.10260 


-37.51 


4.88 


Subaru 


927.06848 


-37.76 


3.84 


Keck 


927.96558 


-28.27 


3.90 


Keck 


931.03690 


-12.10 


3.83 


Keck 


931.94061 


-46.59 


3.92 


Keck 


932.03590 


-42.95 


3.83 


Keck 


932.99985 


3.53 


3.85 


Keck 


933.92455 


66.16 


4.04 


Keck 


934.90368 


49.20 


4.92 


Keck 


934.90692 


49.21 


4.00 


Keck 



for higher Z, with a largely normal distribution. The 
observational error box of ADS 16402B from the SME 
analysis (Fig. |3Jl overlaps with the ZAMS for any Z, and 
naturally, the star cannot lie below the ZAMS for the 
accepted abundance, so this constrains the possible so- 
lution from the lower bounds of log 5 and logL. The 
distribution of possible stellar radii is asymmetric due to 
these constraints and this leads later to asymmetric error 
bars of the planetary radius. 

Our best fit theoretical isochrone is for Z=0.025, and 
roughly bisects the line connecting the centers of the 
SME error boxes. This gives masses of 1.12 ± 0.09 
Mq and 1.16 ± 0.11 Mq and radii of 1.15±|^ J?i?0 and 
1.23±g:};^i?o, for ADS 16402B and ADS 16402A respec- 
tively, and an age of 3.6 Gyrs for the system (see Table 
1), and we use these in the further analysis. 

The absolute / magnitude of ADS 16402B is M/ = 
3.74±0.3 mag (with bolometric correction of -f 0.55 mag). 
From its apparent / magnitude of 9.56 ± 0.06mag (see 
above), and by neglecting reddening (being a near-by 
source at 6 = 20°) we derive a distance of 146±|fpc. The 
same calculation for ADS 16402A yields ISl^^pc. We 
estimate the distance of the ADS 16402AB system as the 
average of the two: 139 ± 20pc. The 11'.'2 separation of 
the two components then implies a projected separation 
of 1550 ± 250 AU. 

All derived stellar parameters, along with absolute and 
apparent / magnitudes, magnitude difference, mean dis- 
tance, and projected separation of the stellar components 
of ADS 16402 are given in Table [2 

4. SPECTROSCOPIC ORBIT OF ADS 16402B 

The Subaru HDS Doppler observations were obtained 
through an Iodine absorption cell (in front of the spec- 
trometer slit for this observing run) to provide a fidu- 
cial wavelength scale and to model th e instrumental 
point spread function fKam be et al.l200^ . We measured 
the star's radial velocity from each high-resolution spec- 
trum, using a synthetic stellar template that was cre- 
ated by modifying a close-match spectrum from a li- 
brary of high resolution t emplate spectra obtained at 
Keck (iJohnson et alJl2006j) . Doppler velocities at Keck 
are also modeled from observations taken through an Io- 
dine cell, however, a standard template observation (ob- 
tained wit hout the Iodine cell ) was used for the reference 
spectrum ijButler et al.lll996|l . Exposure times at both 



Subaru and Keck were typically 500 s, yielding an S/N 
of about 150. 

Emission in the core of the Ca H&K lines was mea- 
sured in each of the Keck spectra to assess the chromo- 
spheric activity of ADS 16402A and B. With five spec- 
tra for ADS 16402A, we measure an average value for 
Shk = 0.16 corresponding to logi?^^ = —4.923 for the 
A component and an average value of nine spectra yields 
Shk = 0.14 and logi?^^ = —5.03 for the B compo- 
nent. These values indicate low chromospheric activity. 
Our estimate for intrins ic astrop hysical velocity jitter is 
3.7m s-i for both stars l|Wrightl f2005). 

Radial velocity measurements from Subaru and Keck 
were used together to derive an orbital fit. The reference 
velocities were those from Keck, and an offset (A^,) be- 
tween Subaru and Keck was also fitted for in the orbital 
solution. Velocity measurements are given in Table [21 
along with observation times and errors. The velocity 
offset of A„ = 14.5m s~^ has been already applied to 
the Subaru data shown in the table, and the expected 
velocity jitter of 3.7m s^^ has already been added in 
quadrature to the uncertainties listed in the table. The 
photometric period of 4.46529 days, and the mid-transit 
time from the photometric ephemerides were fixed in this 
model. We also assumed zero eccentricity based on the 
theoretical expectation of circular orbits for periods as 
short as this. In summary, the fitted parameters were the 
7 velocity of the Keck data points, the K semi-amplitude, 
and the A„ velocity difference between Keck and Sub- 
aru. To determine the uncertainties, 100 Monte Carlo 
fitting trials were run. In each trial, the best fit Kep- 
lerian model was subtracted from the velocities and the 
residuals were shuffied and added back with replacement 
to the theoretical velocities. From these trials, we find 
a velocity semi-amplitude, K = 60.3 ± 2.1ms~^, an rms 
of 5.11ms~^ and a reduced = 1.54. The final fit is 
over-plotted on the RV data in the upper panel of Fig. 21 
The lower panel shows the residuals. 

Adopting a stellar mass of 1.12 ± 0.09Mq, we derive 
Orel = 0.0551 ± 0.0015 AU and msmi = 0.53 ± 0.04Mj. 
Using an inclination ip = 85?9 ± 0?8 (§0 below), this 
yields a planetary mass Mp — 0.53±0.04Afj. The orbital 
elements and uncertainties are summarized in Table O 

There is a small but non-negligible systematic trend 
in the residuals of the orbital fit (lower panel of Fig. 
We note that the fit is improved if the eccentricity and 
oj are allowed to float, and we get e = 0.09 ± 0.02, 
and OJ — 80?7 ± 7?6 with rms of 3.2m s~^ and reduced 
= 0.53. The non-zero eccentricity is only suggestive, 
however, because of the i) small number of data-points, 
and ii) the unrealistically small reduced x^. Nevertheless, 
implications of this finding are potentially very impor- 
tant, and we elaborate further on them in the discussion 

(§I3- 

5. EXCLUDING BLEND SCENARIOS 

We investigated the possibility that the measured 
radial- velocity variations are due not to a planetary com- 
panion but rather to distortions of spectral line profiles 
arising from contamination of the spectrum by the pres- 
ence of an unresolved binary companion to ADS 16402B, 
with an orbital period of ~4.4653 days and an ampli- 
tude of several tens of kms^^. However, this would give 
rise to time-varying asymmetry in the spectral line bi- 
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Fig. 4. — Radial Velocities from Subaru and Keck, folded with 
the period of 4.46529 days. Phase corresponds to the predicted 
time of mid-transit from the orbital ephemeris given in text. Open 
circles represent Subaru data (after the 14.5m s~^ shift described 
in the text, whereas filled circles show Keck measurements. 



TABLE 3 
Parameters of the HAT-P-1b 
planetary system. 



Parameter 



Value 



Period (d) 
T„,id (HJD) 
uj (deg) 
ecc 

Ki (ms-l) 
Ircl (AU) 
m sin i(Mj) 
jp{deg) 
A/p(Mj) 

RpjRj) 



4.46529 ± 0.00009=^ 
2453984.397 ± 0.009=^ 

Qa 
Qa 

60.3 ±2.1 
0.0551 ±0.0015 
0.53 ±0.04 

85.9 ±0.8 
0.53 ±0.04 

1 Sfi+O.ll 



Fixed in the orbital fit. 



sectors (see, e.g.. ISantos et 7ni2nnl iTorres et a].ll2nnl. 

We have searched for such Une bisector variations in our 
Subaru spectra and find no evidence for them. Compar- 
ison between the relative radial velocities and bisector 
spans of ADS 16402A and B is shown in Fig. The 
results are based on the Subaru spectra only, because 
they were obtained on the same nights for both stars. 
Line bisectors for each spectrum were computed from 
the cross-correlation function averaged over all orders, 
which is representative of the average spectral line pro- 
file of the star. The cross-correlations were performed 
against a synthetic spectrum matching the properties of 
the stars. The bisector spans were then computed as 
the velocity difference between poi nts chosen near th e 
top and bottom of the bisectors (see ITorres et a,1.ll2nnl . 
ADS 16402A shows no variation in either the velocities 
or the bisectors spans. ADS 16402B, on the other hand, 
shows obvious velocity variations but does not show the 
concomitant change in the bisector spans that would be 
expected if it were the result of a blend with an eclips- 
ing binary. If the light and velocity variations were the 
result of a blend with an eclipsing binary, the bisector 
spans would be expected to vary with a similar ampli- 
tude as the velocities (see, e.g.. lQueloz et al.l200H) . That 
scenario is therefore ruled out by these observations. 

We have also modeled th e light curve as a blend fol- 
lowing ^b^^^3 l|2004b|) . subject to all available ob- 
servational constraints. Specifically, we calculated model 
light curves for a three-star system consisting of a slightly 
evolved early G star primary plus a later-type secondary 
and a tertiary which together comprise an eclipsing bi- 



FlG. 5. — Comparison between the relative radial velocities 
and bisector spans of ADS 16402A (left panel) and ADS 16402B 
(right panel), based on the Subaru spectra. ADS 16402A shows no 
variation in either the velocities or the bisectors spans (left panels), 
while ADS 16402B shows obvious velocity variations but does not 
show the variation in the bisector spans that would be expected 
if the velocity signal were the result of a blend with an eclipsing 
binary. 

nary, with the requirement that the brightness of the 
secondary be no more than '^10% that of the primary 
(to have escaped detection in the DS spectra) and that 
the light curve of the eclipsing binary, attenuated by the 
much brighter primary, should be indistinguishable from 
the observed light curve in the FLWO 1.2 m data. How- 
ever, we have found it impossible to construct an accept- 
able light curve in this fashion. 

Based on the above, we conclude that ADS 16402B is 
orbited by a transiting object that gives rise to the ob- 
served photometric and radial velocity variations. From 
the 60m s~^ amplitude of the radial velocity orbit (§0, 
we conclude that the object must have mass about 0.53 
Mj and hence must be an extrasolar planet. 



6. 



THE TRANSITING PLANET HAT-P-lb AND ITS 
PROPERTIES 



In order to determine more precisely the physical char- 
acteristics of HAT-P-lb and its orbital parameters, we 
have carried out model fitting of the different sets of pho- 
tometric data shown in the four panels of Fig.|21 A simple 
goodness-of-fit was determined against t he predicted flux 
values evaluated with the algorithm of iMandel fc Agol 
(2002'). We used the values for the stellar radius and 
mass (and their uncertainties) determined above, and co- 
efficients for stellar limb darkening for the Sloan z-filter 
by Claret (2004). First we fitted the original discovery 
data from HATNet with the appropriate correction for 
the flux of the binary star (unresolved on the HATNet 
images): the Rp = 1.42i?j and ip = 85?4, have large 
uncertainties, as anticipated, and might be subject to 
large systematic errors. We then fitted to the recent 
Konkoly Schmidt observations of a complete HAT-P-lb 
transit in /-band: the values are Rp = 1.31±0.07i?j and 
ip = 86?0±1?0, and a precise estimate of transit duration 
of d = 0.11671 days and center of transit = 3984.3967. 
The latter value allows us to estimate accurately the 
of the transit we observed just an orbital period earlier 
(Aug 31) at FLWO, where the beginning of the ingress is 
missing (see Fig. |21 third panel from the top) . Combin- 
ing both z-band transit curves from FLWO with a known 
Tc allows an independent fit: we get Rp = 1.36±0.05i?j 
and ip = 85?9 ± 0?8. The fit is better because the qual- 
ity and amount of photometric data is higher, and the 
limb darkening in z-band is less. The fits to all data 
sets give us the same planet parameters within the error; 
henceforth we use our best fit with the FLWO z-band 
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transits. 

The result for the planetary radius is Rp = 
1.36±Q jgi?j, where the above statistical errors from the 
fit and systematic errors from the uncertainties of the 
stellar radius and mass (ztQ-gy), have been added in 
quadrature. Due to the asymmetric error bars on the 
stellar radius (see above) we have asymmetric error bars 
on Rp too. The orbital inclination is ip = 85?9 ± 0?8. 
The radius, orbital inclination, and derived mass of the 
planet are included in Table 13 

An even better light curve with e.g. the Hubble Space 
Telescope, and a more detailed analysis of the parameters 
of the stellar pair could reduce that error considerably. 
The results of this on-going work will be reported in a 
forthcoming paper, as soon as the multi-band observa- 
tions enable us to refine the binary system parameters. 

7. DISCUSSION 

Fig. shows the location of HAT-P-lb in a plot of 
planetary radius versus planetary mass for the now 11 
known transiting exoplanets, plus Jupiter and Saturn. 
It is immediately seen that HAT-P-lb is in an extreme 
position. It has a radius apparently as large as or larger 
than the largest known exoplanet, HD 209458b, and at 
the same time it has a mass that is significantly smaller. 
Thus, it has a lower mean density than any other known 
planet. Until now, HD 209458b was the only signifi cantly 
anomalous transiting planet. Thus. iGuillotl 1)2005(1 noted 
that all of the transiting planets known at that time, with 
the exception of HD 209458b, fall within a mass-radius 
relation derived from modeling the interiors of strongly 
irradiated planets. All of the transiting planets detected 
subsequently, except HAT-P-lb, also fit the same rela- 
tion. Among them, OG LE-TR-lOb had been thought to 
be similarly anomalous fKonacki et al."2005'), but recent 
evidence (Holnian ct al. 2005; Santos et al. 200 6. M. Hol- 
man, private communication) points to a smaller plan- 
etary radius an d normal density. Simil arly, the XO-lb 
discovery paper l|McCullough et al.l2006|) derives a larger 
radius than the value later refined by l|Holman et al.l 

iooe). 

HD 209458b has a radius some 10 to 20 percent larger 
than expected, given its age, mass, and expected temper- 
ature including effects of irradiation. Possible solutions 
include: (a) a fraction of a percent of the stellar flux irra- 
diating HD 209458b could be transported to deep layers 
of order 1 00 bar, and thereby keep the planet's radius 
puffed up l)Showman fc Guillotl2002D : (b) there is a small 
eccentricity of its orbit, forced by the presence of a so-far 
undetected planetary companion and leading to internal 
tidal heating (Bodcnheimcr. Laughlin. & Lin 2003), and 
(c) HD 209458b might be in a Cassini State — a dynami- 
cal state where its rotation axis lies close to the plane of 
its orbit l|Winn fc Holmanl l2b05). leading to large obliq- 
uity tides that cause heating at significant depth leading 
to the planet's large radius. A natural question about 
the first of these mechanisms is why all close-in "hot 
Jupiters" are not similarly bloated by that mechanism. 
A problem with the last two mechanisms is that both 
are ad hoc] they require unusual circumstances for which 
there is no independent evidence. 

With the discovery of another planet, HAT-P-lb, 
which has a radius similar to and probably even larger 
than HD 209458b, and an even smaller density, mecha- 




/ /■ I''l49626b 

/ ,/•■ J- ' [p]: g cm" 

0.6 - ; 



0.5 1 1.5 2 

Mpl [Mjupl 

Fig. 6. — Mass-radius diagram depicting known tran- 
siting exoplanets, plus Saturn and Jupiter (for comparison). 
Sources for the mass and radius values (in this order) are 
listed after the name of the planet; HAT-P-lb: th is work , 
HD 209458b: iLaughlin et al. (200?), "Knutson et al.i f200^ 
HD 18 9733b: IBakos et alT l2006'). OGLE-TR-lllb: Pont ct al] 
J^Ol), OGLE-TR-lOb: (Santos ct al. (2006), OGLE-TR-132b: 
"Moutou ct al. 12004), OGLE-TR-56b: Torres ct al. 1 2004a), 
OGLE-TR-113b: Bouchy et al. (2004); Konacki et al. (2004), 
TrES-1: Sozzctti ct al. |200|i), HD 149026: Sato ct al. (2005), 
l^harbonncau ct al. 1 2006 a), X O-lb: iHolman ct al. (2006), 
TrES-2: lO'Donovan et all i2?OT) . 

nisms requiring unusual circumstances seem less likely. 
Moreover, for HD 209458b, ex planation (b) appears 
ruled out because Doppler data l|Laughlin et al.l 12005131 
iWim^^y [2QQ ^ and transit timing variations data 
l|Miller-Ricci et al."200?) have excluded companion plan- 
ets capable of enacting it. 

In the case of HAT-P-lb, it remains possible that the 
system has a second planet capable of forcing its eccen- 
tricity and leading to internal tidal heating. Our current 
orbital solution (§0)1 suggests e = 0.09 ± 0.02. If con- 
firmed, such eccentricity would be a strong indication 
for the presence of a third body (second planet), because 
otherwise tidal dissipation would have circularized HAT- 
P-lb's orbit in ~ 225( j^) Myr. The dissipation factor, 
Q, is poorly known for "hot Jupiters", but is expected 
to lie between the values for Jupiter and solar-like stars 
(10^ to 10*^). With forced orbital eccentricity, e = 0.09, 
HAT-P-lb will be subject to internal tidal heating due 
to an energy dissipation rate 

dE _ me^n fRpV' GMj 
dt ~ 4 Q \ a J a 

where G is the gravitational constant, M^, is the mass 
of the star, a is the semimajor axis, n is the mean 
motion, and Q is the specific dissipation factor. With 
e = 0.09, the rate of energy dissipation in HAT-P-lb 
would be 3.0 x IQ^^erg s~^, corresponding to energy dis- 
sipation per unit Jupiter mass of 5.7 x lO^^erg s~^M^^ . 
This rate is almost identical to the dissipation rate 
per unit Jupiter mass of 5.8 x lO^^erg s^^MJ^ re- 
quired to achieve a 20% increase in i?p for HD 209458b 
(no core model) and explain its a nomalous low den- 
sity l|Bodenheimer. Laughlin. &: Linll2003D . Therefore we 
conclude that if the orbital eccentricity is confirmed, the 
anomalous Rp and low density of HAT-P-lb could be 
easily explained by tidal heating. The perturbing sec- 
ond planet would be detectable by both transit timing 
variations and Doppler velocities. 

For now, however, all that can be said is that bloated 
planets are not that unusual, and, since a second perturb- 
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ing planet is lacking in HD 209458b's case, perhaps some 
other explanation should be sought that can explain such 
planets more generally. 
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